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Abstract We study shear-free spherically symmetric relativistic models with heat 
flow. Our analysis is based on Lie's theory of extended groups applied to the 
governing field equations, fn particular, we generate a five-parameter family of 
transformations which enables us to map existing solutions to new solutions. All 
known solutions of Einstein equations with heat flow can therefore produce infinite 
families of new solutions. In addition, we provide two new classes of solutions 
utilising the Lie infinitesimal generators. These solutions generate an infinite class 
of solutions given any one of the two unknown metric functions. 
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1 Introduction 

In this paper, we consider spherically symmetric radiating spacetimes with van- 
ishing shear which are important in relativistic astrophysics, radiating stars and 
cosmology. The assumption that the shear vanishes in a spherically symmetric 
spacetime, in the presence of nonvanishing heat flux, is often made to describe the 
dynamics of cosmological models. The importance of relativistic heat conducting 
fluids in modeling inhomogeneous processes, such as galaxy formation and evolu- 
tion of perturbations, has been pointed out by Krasinski [1]. Some of the early 
exact solutions in the presence of heat flow were given by Bergmann [2] , Maiti [3] 
and Modak [4] . Deng [5] provided a general method of generating solutions to the 
Einstein field equations which contains most previously known exact solutions. 
Heat conducting exact solutions are necessary to generate temperature profiles 
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in dissipative processes by integrating the heat transport equation as shown by 
Triginer and Pavon [6]. Bulk viscosity with heat flow affects the dynamics of in- 
homogeneous cosmological models as shown by Deng and Mannheim [7] . Recently 
Banerjee and Chatterjee [8] and Banerjee et al [9] have investigated heat conduct- 
ing fluids in higher dimensional cosmological models when considering spherical 
collapse, the appearance of singularities and the formation of horizons. The role of 
heat flow in gravitational dynamics and perturbations in the framework of brane 
world cosmological models has been highlighted by Davidson and Gurwich [10] 
and Maartens and Koyama [11]. 

The presence of heat flux is necessary for a proper and complete description 
of radiating relativistic stars. The result of Santos et al [12], in one of the first 
complete relativistic radiating models, indicates that the interior spacetime should 
contain a nonzero heat flux so that the matching at the boundary to the exterior 
Vaidya spacetime is possible. Models containing heat flow in astrophysics have been 
applied to problems in the gravitational collapse, black hole physics, formation 
of singularities and particle production at the stellar surface in four and higher 
dimensions. The study of Chang et al [13] showed that the process of gravitational 
collapse of a spherical star with heat flow may serve as a possible energy mechanism 
for gamma-ray bursts. 

Herrera et al [14], Maharaj and Govender [15], and Misthry et al [16] have 
shown that relativistic radiating stars are useful in the investigation of the cos- 
mic censorship hypothesis and in describing collapse with vanishing tidal forces. 
Wagh et al [17] presented solutions to the Einstein field equations for a shear- free 
spherically symmetric spacetime, with radial heat flux by choosing a barotropic 
equation of state. For particles in geodesic motion a general analytic treatment 
is possible and solutions are obtainable in terms of elementary and special func- 
tions as demonstrated by Thirukkarash and Maharaj [18]. Herrera et al [19] found 
analytical solutions to the field equations for radiating collapsing spheres in the 
diffusion approximation. These authors demonstrated that the thermal evolution 
of the collapsing sphere which can be modeled in causal thermodynamics requires 
heat flow. Note that stellar models with shear are difficult to analyse but partic- 
ular exact solutions have been found by Naidu et al [20] and Rajah and Maharaj 
[21]- 

Shear-free fluids are also essential in modeling inhomogeneous cosmological 
processes. Krasinski [1] points out the need for radiating models in the formation 
of structure, evolution of voids, the study of singularities and cosmic censorship. 
Heat conducting fluids are important in cosmological models in higher dimensions 
and permits collapse without the appearance of an event horizon; this aspect has 
been studied by Banerjee and Chatterjee [8]. In brane world models the presence 
of heat flow allows for more general behaviour than in standard general relativity, 
the analogue of the Oppenheimer-Snyder model of a collapsing dust permits a 
radiating brane which was proved by Govender and Dadhich [22]. 

In this paper we intend to analyse the pivotal equation previously studied by 
Deng [5] . He developed a general (though ad hoc) method to generate solutions and 
obtained a new class of solutions which included various known special cases (see 
Sect. 2.). We adopt a systematic approach (using Lie theory) to generalise known 
solutions and generate new solutions of the same equation. The basic features of 
Lie symmetry analysis are outlined in Sect. 3. In Sect. 4, we extend the Deng [5] 
known solutions to find new solutions to the fundamental equation utilising Lie 
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theory. In Sect. 5, we systematically study other group invariant solutions admitted 
by the fundamental equation. For most of the symmetries, we obtain either an 
implicit solution or we can reduce the governing equations to a Riccati equation 
which is difficult to solve in general (though particular solutions can always be 
found) . There are two cases in which we find new exact solutions regardless of the 
complexity of the generating function chosen. Some brief concluding remarks are 
made in Sect. 6. 



2 Radiating Model 



Due to the requirements of spherical symmetry and the shear-free condition, the 
line element can be written as 

1 



ds 



-D 2 dt 2 + |dr 2 + r 2 {d6 2 + sin 2 ( 



(1) 



where D and V are functions of t and r. In the study of solutions of the Ein- 
stein equation with heat flux, Deng [5] studied a shear-free spherically symmetric 
cosmological model where he considered the energy-momentum tensor as 



T^ v = {p + p)U^U v +pg^v + q^Uu + qvUfi 



(2) 



where Up, p, p are the four- velocity of the fluid, energy density and pressure, 
respectively, and is the heat flux. The Einstein field equations are given by 
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which is the condition of pressure isotropy with u = r 2 . Glass [23] and Bergmann 
[2], also discovered that in the comoving system, Einstein field equations generate 
the pressure isotropy condition given by the equation (4), which is the master 
equation for the system (3a)- (3d). 

A number of authors have obtained various solutions to (4), among which is 
the conformally flat class 

c{t)u + d{t) 



Equations (36) and (3c) together imply 

VD UU + 2D u Vu - DVuu 



D a(t)u + b(t) 
V = a(t)u + b(t) 



(5a) 
(5b) 
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where a, b, c, and d are arbitrary functions of t. Sanyal [24] and Modak [4] obtained 
this class along with other solutions, while Bergmann [2] and Maiti [3] obtained 
special cases of the class. 

A method of generating more general solutions to the master equation (4) 
has been developed by Deng [5] who found solutions when simple forms of V or 
D are chosen. In finding solutions, Deng [5] considered the master equation as an 
ordinary differential equation with respect to u. He treated (4) as a linear equation 
of V if D is a known function and vice versa. His approach was as follows: 

— Choose a simple function D = D\ and find the most general solution V = V\ 
of (4). 

— Take V = V\ and find the most general solution D = D2 obeying (4). 

— Take D = D2 and find the most general solution V = V2. 

This procedure can be continued indefinitely. By alternating between D and V, this 
process can go on forever generating infinite series of solutions expressed in terms 
of integrals. This is a powerful method as all known solutions can be generated 
using this algorithm 

In this paper, we show the Deng [5] general method of generating solutions, 
that gives a general class of solutions which include (5a)-(5b) as special cases, 
may be extended by a simple invariant transformations. In addition, we reduce 
the order of (4) via Lie analysis to obtain new solutions not obtainable via the 
Deng approach. 



3 Lie analysis 

The basic feature of Lie analysis for a system of ordinary differential equations in 
two dependent variables requires the determination of the one-parameter (e) Lie 
group of transformations 

u = f{u,V,D,e) (6a) 
V = g(u,V,D,e) (6b) 
D = h(u, V, D, e) (6c) 

that leaves the solution set of the system invariant. It is difficult to calculate these 
transformations directly. As a result, one tends to look for the infinitesimal form 
of the transformations, viz. 

u = u + e£(u,V,D)+0(e 2 ) (7a) 
V = V + er,{u,V,D)+0(e 2 ) (7b) 
D = D + sC(u,V,D) + 0(e 2 ). (7c) 

This form of the transformations can be obtained once we obtain their "generator" 



G = ^ +11 W +(: dD 



(8) 
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(also called a symmetry of the differential equation) which is a set of vector fields. 
Having found the symmetries, we can regain the finite (global) form of the trans- 
formation (6a)- (6c), by solving 

^=t{u,V,D) (9a) 
^=n(u,V,D) (9b) 
^ = t(u,V,D) (9c) 



subject to 

u| e =o=«, V\ e =o = V, D\ £=0 = D. (10) 

(The full details can be found in a number of excellent texts such as Bluman and 
Kumei [25] and Olver [26]). 

The determination of the generators is a straight forward process, greatly aided 
by computer algebra packages (see the treatments of Dimas and Tsoubelis [27] 
and Cheviakov [28]). We have found the package PROGRAM LIE by Head [29] to be 
the most useful in practice. It is quite remarkable how accomplished such an old 
package is; it often yields results when its modern counterparts fail. 

Utilising PROGRAM LIE, we can demonstrate that (4) admits the following Lie 
point symmetries/vector fields: 



Gi = — 
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(11a) 


G 2 = u— 
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At this stage, it is usual to use these symmetries to reduce the order of the 
equation in the hope of finding solutions. Before we proceed with this approach, 
we show how the finite form of the transformations generated by these symmetries 
can generate new solutions from known solutions. 



4 Extending known solutions 

Having found the symmetries of (4) we know that they generate transformations 
of the form (6a)- (6c) that leaves (4) invariant. 

We illustrate the approach using the infinitesmal generator G\ . First we observe 
that 



C = l, tj = 0, C = 



(12) 
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We solve equations (9a)-(9c), subject to (10), to obtain 

u = u + a\ (13a) 
V = V (13b) 
D = D (13c) 

This means that (13a)-20cmaps equation (4) to the form 

VDuu + 2DuVu - DVuu = 0. (14) 

As a result, any existing solution to equation (4) can be transformed to a solution 
of (14) (and so a solution of (4) itself) by (13a)- 20c. Note that, usually, a\ is an 
arbitrary constant. However, since V and D depend on u and t we take oi to be 
an arbitrary function of time, a\ = ai(t). 

If we now take each of the remaining symmetries successively, we obtain the 
general transformation as follows: 

(15a) 
(15b) 
(15c) 
(15d) 

(15e) 

The combination of all the transformation of symmetries, therefore leads to the 
general relationship: 



Gi : 


: u = 


a\ + u, 


D 


= D, 


V = 


V 


G 2 


: u = 


e a2 { ai +u), 


D 


= D, 


V = 


V 


G 3 


: u = 


e a2 { ai +u), 


D 


= e a3 D, 


V = 


V 


d 


: u = 


e a2 { ai +u), 


D 


= e a3 D, 


V = 


e ai V 


G 5 : 


; u = 


e Q2 (ai +u) 


D 


= e a3 D, 


V = 


e ai V 


1 - a 5 e a 2(ai + it) ' 


1 — a 5 e a2 (ai + u) 



oi + u) 



(16a) 

D = e a3 D (16b) 
e ai V 

V = —r 7 (16c) 

1 - a 5 e a 2 (oi + u) ' 

where the a, are all arbitrary function of time. 

Thus any known solution of equation (4) can be transformed to a new solution 
of equation (4) via (16a)-(16c). For example, the particular Deng [5] solution 

D = 1, V = a(t)u + /3(t) (17) 

is transformed to the new solution 

e a2 (ai +«) 

u = — -. — r (18a) 

1 - a 5 e a 2 (oi + u) K ' 

v = pmt+m (18b) 

D = e a3 (18c) 

All the solutions in the Deng [5] class, the conformally fiat models (5a)-(5b), the 
result listed in Krasinski [1] and Stephani et al [30] can be extended by (16a)-(16c) 
to produce new solutions of (4). Also note that all the new results that we derive 
in the next section can be similarly extended via (16a)-(16c). 
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5 New Solutions via Lie symmetries 

The usual use of symmetries of ordinary differential equations is to reduce the 
order of a differential equation. For symmetries (lla)-(lle) we obtain either an 
implicit solution of (4) or we can reduce the governing equations to a Riccati 
equation. Both these results are no improvement to that of Deng [5] , i. e. we need 
to choose simple forms for one of the functions (either D or V) in order to solve 
for the other. However there are two cases in which we can find new solutions 
regardless of the complexity of the function chosen. 



5.1 The choice D = D(V) 

An obvious relation to consider is when one dependent variable in (4) is a function 
of the other. In general such an approach usually results in a more complicated 
equation. However, using the Lie symmetry Gi (which gives the same result as 
G2), we can make significant progress. For our purposes we use the partial set of 
invariants of 

Gl = I < 19 > 

given by 

p = V (20a) 
q(p) = V u (20b) 
r(p) = D (20c) 

This transformation reduces equation (4) to 

q'{p) [r{p) - pr'(p)] = q(p) [pr"(p) + 2r'(p)] (21) 

which can be integrated to give 

q = qoe J r-r'p dp (22) 

Substituting for the metric functions via (20a)- (20c), we can integrate one more 
time to give the solution 



/ 



dV = qou + uq (23) 



where go an d wo are arbitrary functions of time. This means that, given any func- 
tion V depending on D, we can work out V explicitly from (23). Such a relationship 
between V and D has not been found previously Note that since (4) is linear, once 
we obtain V via (23) we can use it to obtain the general solution of (4) using stan- 
dard techniques for solving linear equations. (In the case of G2, we obtain exactly 
the same solution as (23).) 

We illustrate this method with simple examples. Using the particular Deng [5] 
solution D = 1, we evaluate (23) to obtain 

V = q (t)u + u (24) 



which is exactly what Deng [5] obtained. 
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If we take D = V 2 , then equation (23) is reduced to 

j V 6 dV = q u + u (25) 

and hence 

Vi = {q u + u ) 1/7 (26a) 
D = {q u + u ) 2/7 (26b) 

which is a solution of (4). Having found V from (23) we can generate a new solution 
as follows: the second independent solution of (4) has the form 

V 2 = y(qou + u ) 1/7 (27) 
with y being unknown. Then (4) becomes 



2q 



y 



7 (q u + mq) 



y' = (28) 



with solution 



Hence 



y=^(q u + u ) 9/7 + C 2 (29) 
qo 

V = ^ (q u + u ) 10/7 + C 2 (qou + u ) 1/7 (30) 

qo 

where Ci and C 2 are arbitrary functions of time, is the general solution to (4) 
when D = (qou + uq) 2 ^ 7 . 



5.2 The choice W = % 

We also consider the ratio of V to D (and later D to V) to generate a new solution. 
Incidentally both ratios arise as a result of a combination of the generators G3 and 
G4 given by 

Gs + G 4 = D^ + V^ (31) 
This symmetry combination gives rise to the invariant 

W = J (32) 
Then equation (4) is transformed by (32) to the form 

-2WD 2 u + D 2 Wuu = (33) 

with solution 

D = C 1 ( i )exp(|±^dn) (34) 

Given any function W = W(u) we can integrate the right hand side and find a 
form for D. If we take W = a(t)u + b(t), then (34) gives 

D = Ci(t) (35) 
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and 

V = C 1 (t)(a(t)u+b(t)) (36) 

which corresponds to a Deng [5] solution. 

Alternatively, we could substitute the inverse of (32), i.e. 

W = y (37) 

into (4) and obtain 

2D I2 W 2 - 2D 2 W' 2 + D 2 WW" = (38) 




W' 2 - hww \ 

du (39) 



with solution 



Again, given any function W = W{u) we can integrate the right hand side and 
find a form for D. If we take W = a(t)u + b(t) as before, in (39), we find that 

D 1 =C 2 {t){a(t)u+b{t)) (40) 

and 

Vt=C 2 (t) (41) 
which again is essentially a solution of Deng [5]. However, we also have 

D2 = (42) 

and 

V 2 = ^ (43) 

satisfies (4), thus obtaining two different solutions from the same seed function. 

Observe that (34) and (39) will contain all solutions of Deng [5] for appro- 
priately chosen seed functions W or W which are ratios of the metric functions. 
While Deng's approach requires simply chosen forms of either D or V in order to 
integrate (4), we have no such requirement. We are always able to reduce (4) to 
the quadratures (34) or (39) regardless of the complexity of the seed functions, a 
result not obtainable via Deng's approach. 



6 Conclusion 

We have been able to extend Deng's solutions [5] of the Einstein field equations 
governing shear-free heat conducting fluids in general relativity. This was accom- 
plished by using simple transformations based on the invariance properties of the 
equation under study viz. (4). In addition, motivated by the invariants of the 
symmetries admitted by (4) we were able to reduce (4) to quadratures for any 
given seed function. This leads to three new classes of solutions for infinite fam- 
ilies of functional forms involving D(V), W and W. This is an improvement on 
the approach of Deng [5] who required 'simple' functional forms for D or V to 
be chosen before (4) could be solved. This again promotes the use of Lie's the- 
ory of extended groups to analyse the Einstein field equations arising in different 
applications/models in general relativity. 
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